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ABSTRACT 

The composition and energy spectra of charge 
groups (C - 0), (Ne - S), and (Z ^ 17) above 500 
GeV/nucleon from the experiments of JACEE series 
balloon borne emulsion chambers are reported. 

1. Introduction . 

Studies of cosmic ray elemental composition at higher energies 

provide information on propagation through interstellar space, accelerat- 

ion mechanisms, and their sources. One of the present interests is the 
elemental composition at energies above 100 GeV/nucleon. Statistically 

sufficient data in this energy region can be decisive in judgement of 

propagation models from the ratios of SECONDARY/PRIMARY and source spectra 

(acceleration mechanism), as well as speculative contributions of differ- 
ent sources from the ratios of PRIMARY/PRIMARY. At much higher energies, 

i.e., around 10 15 eV, data from direct observation will give hints on the 

"knee” problem, as to whether they favor an escape effect possibly 

governed by magnetic rigidity above 10 1 6 * eV. 

The JACEE balloon flight experiments continue to measure composition 

and energy spectra of cosmic rays (Z = 1 to 26) directly at energies 10 12 * 
to 10 15 eV using large area thin emulsion chambers [1 - 5]. The previous 
JACEE results indicated no significant change of spectral indices up to 

500 TeV and 50 TeV/nucleon for proton and helium spectra, respectively [4]. 

Flux values of each group (C - 0), (Ne - S), and (Fe) at least up to 

lO 14 eV also indicated no significant evidence for heavy nuclei dominance 
within the limited statistics [5]. In this paper, the results for the 



49 


OG 4.1-13 

energy spectra of charge groups (C - 0), (Ne - S), and (Z ^ 17) above 
500 GeV/nucleon are updated. 

2.. Apparatus and Experimental Procedure . 

” Details of the apparatus and techniques of JACEE emulsion chambers 
have been reported in refs. [1 - 5]. Each detector, comprised more than 
350 layers of materials, is functionally divided into three sections: (1) 
the pri mary charge determination module at the top, (2) the target module 
where nuclear interactions occur preferdntialy (3) the calorimeter at the 
bottom to measure the energies of released gamma rays from interactions. 

Events were detected by visual scanning of X-ray films for dark spots 
produced by electromagnetic cascades in the calorimeter under the same 
criterion of darkness (D) in an area of 200 x 200 ym 2 . For these events, 
primary charge Z (AZ < 1) above the interaction and the total energy of 
released gamma rays zEy (AZEy < 25 %) were measured [1 - 5], 

The maximum value of darkness (D m j in the cascade development is a 
function of not only zEy but also the vertex height (H) from the top of 
the calorimeter, incident angle (e), and interaction characterstics 
(transverse momentum and multiplicity). Among these parameters, the main 
factors are z£y and H. For a fixed zEy value. Dm distributes around some 
mean value. The width of this distribution g(Dm)dDm depends on the 
primary charge (mass) and the chamber design. g(Dm) has been calculated 
by a Monte Carlo method using the chamber design adopted here, interaction 
characteristics (CKP and superposition models), and characteristic curves 
of the response of X-ray films. Then, the detection efficiency P(Z, ZEy) 
under the threshold. D^' 1 is calculated by, 

P(Z, zEy) = / nth g(Gm)dDrn / /n g(Dm)dDm . 

>u m u m 

In case of calorimeter jets, i.e., H = 0, P(zEy) rapidly reaches 100 % at 
zEy = 1 TeV for all elements under the threshold level D# 1 =0.1 which is 
a practical threshold, because of no fluctuation of H. On the other hand, 
in the case of target jets, the practical threshold becomes about 0.15 
and the width of g(Dm) becomes larger due to a fluctuation of H values. 
This gives higher values of zEy to obtain P(zEy) = 100 %. The critical 
zEy values have been estimated by the above method for different nuclear 
species. For light elements such as proton and helium, the critical zEy 
was about 1.5 TeV which has been achieved in an analysis of selected 
JACEE chambers. Under the present standard event selection (D^ h = 0.3) 
the detection efficiency P(zEy ^ 10 TeV) is confirmed to be 100 % for all 
elements and for all existing JACEE chambers. 

3. Deconvolution of Primary Spectrum. 

The observed zEy spectrum is a convolution of the primary spectrum 
and the distribution function f ( ky ) where ky (= zEy/Eo) is a partial 
inelasticity into gamma rays; Eo is primary energy. f(ky) has been also 
calculated by the above Monte Carlo method [63, including succesive inte- 
ractions, assuming its energy independence. The ZEy spectrum F(zEy)dzEy 
is given by, 

F(zEy )dzEy = f ky f(ky)dky f £ G(E 0 )dE 0 «(zEy - kyEo) , 

where G(Eo)dEo is primary spectrum. This formula directly shows that the 
ZEy spectrum is uniquely a power law with the same index as the primary 
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spectrum when G(E 0 )dE 0 « E 0 “ 6 dE 0 . The energy conversion factor Cky from 
the zEy spectrum to primary spectrum. is then given by, 

Cky = [ / k ky ^- 1 f ( ky ) dky H 1 /^' 1 ) 

Cky is applied to obtain the same flux between the zEy spectrum and the 
primary spectrum at Eq = zEy/Cky. Values of Cky for each element in case 
of calorimeter and target jets were calculated as follows, 

P He N Mg Ca Fe 

calorimeter 0.240 0.174 0.100 0.082 0.064 0.058 

target 0.222 0.159 0.120 0.117 0.115 0.115 

for chambers mainly used in this work [71. 

The collecting power of the detectors was calculated by a Monte Carlo 
method using both the geometric aperture and the probability of interac- 
tions for each element. Absolute fluxes above chambers were estimated 
considering the resolution function of zEy measurements which is an almost 
constant Gaussian type with a 25 % width in our case. After atmospheric 
correction at the depth of 3 - 5 g/cm 2 , the flux values at the top of 
atmosphere were obtained. 

4. Results and Discussion. 

Figs. 1, 2, and 3 show the results of the integral energy spectra for 
groups (C - 0), (Ne - S), and (Z ^ 17) above 500 GeV/nucl eon under the 
present event selection criterion of D£ h based upon a simulation. In 
Fig. 3, Fe spectrum at energies 20 to 80 GeV/n from the JACEE-3 hybrid 
detector experiment [8] is also shown. In these figures, the low energy 
results from two experiments [9, 10], recently reported, are included for 
comparison. Solid lines are the extrapolation from the flux value at 25 
GeV/n of ref. [9] with the spectral indices 1.7 and 1.5. Data around 4 
TeV/n for (C - 0) and (Ne - S),and those around 1 TeV/n for (Z ^ 17) are 
corrected, by empolying the Monte Carlo calculations for data analysis 
below the sub-threshold zEy = 10 TeV. In these figures, at energies 
higher than sub-threshold, data of each group indicate slightly higher 
flux values (to about one s.d. from the extrapolated spectra with an index 
1.7) with the statistics of only 4, 2, and 1 of 3 (Ca at 78 TeV/n) events 
for respective groups, while (C - 0) data around 1.5 TeV/n with the stati- 
stics of 10 events are in agreement with the low energy data for index 1.7. 
Statistical improvements are obviously desired for our continueing efforts 
at the highest energies above 10 14 eV. 

Can emulsion chamber technique ident'fy a bending of the primary spe- 
ctrum by using of the zEy method under the convolution of ky distribution? 
Monte Carlo studies have shown a positive answer for this question, at 
least for heavy nuclei component. If materials of more than about 3 
interaction length are available, then the distribution function f(ky) 
approaches the Gaussian type with about 35 % width, by virtue of succesive 
interactions. Hence, it is sufficient to detect any narrow peaks or bends 
in the energy spectrum, provided sufficient statistics are gathered. For 
observations above 10 15 eV, however, the energy resolution becomes poorer 
than our current observations, because the shower maximum of electro- 
magnetic cascades in the calorimeter becomes too deep for thin chambers 
currently adopted. Even in this case, emulsion chamber techniques still 
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provide measurements of cosmic ray energy above 10 15 eV without increasing 
the thickness of calorimeter if the e* pair method is used to augment the 
zEy energy determination [11 ]. 
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